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The numerically exact superposition T-matrix method is used to model far-field electromagnetic scattering by two 
types of particulate object. Object 1 is a fixed configuration which consists of N identical spherical particles (with N 
= 200 or 400) quasi-randomly populating a spherical volume V having a median size parameter of 50. Object 2 is a 
true discrete random medium (DRM) comprising the same number N of particles randomly moving throughout V. 
The median particle size parameter is fixed at 4. We show that if Object 1 is illuminated by a quasi-monochromatic 
parallel beam then it generates a typical speckle pattern having no resemblance to the scattering pattern 
generated by Object 2. However, if Object 1 is illuminated by a parallel polychromatic beam with a 10% bandwidth 
then it generates a scattering pattern that is largely devoid of speckles and closely reproduces the quasi- 
monochromatic pattern generated by Object 2. This result serves to illustrate the capacity of the concept of 
electromagnetic scattering by a DRM to encompass fixed quasi-random particulate samples provided that they are 
illuminated by polychromatic light. 
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1. INTRODUCTION 


A discrete random medium (DRM) is defined as a 
morphologically complex object in the form of an imaginary 
volume V populated by a large number N of small particles in 
such a way that the spatial distribution of the particles 
throughout the volume is statistically random and uniform [1- 
5]. Typical examples of DRMs are natural and manmade 
suspensions of particles in gases and liquids. Of course at a 
given moment in time, the spatial distribution of particles in 
any multi-particle group is definite rather than random. 
Therefore, if the group is illuminated by a monochromatic or 
quasi-monochromatic parallel beam of light then statistical 
randomness and spatial uniformity of a DRM can be achieved 
only over a sufficiently long period of time as a result of random 
temporal changes of particle positions. This is precisely what 
occurs naturally in a multi-particle group suspended in a gas or 
a liquid and results in smooth, speckle-free patterns of time- 
averaged electromagnetic scattering. 

The concept of a DRM is also frequently used in application 
to fixed particulate media such as powder surfaces, sheets of 


paper, layers of paint, or biological tissues. This usage is usually 
motivated by the non-detection of speckles in measured 
scattering patterns. It should be recognized, however, that in 
such cases the non-detection of speckles is explained by factors 
other than random movements of particles, for example, by the 
illumination of a fixed multi-particle group by a polychromatic 
beam [6]. 

Yet the smoothness of the scattering patterns measured for 
fixed particulate media often serves as a justification for using 
theoretical modeling techniques developed for true DRMs. 
Among such techniques are the theories of radiative transfer 
and weak localization derived from the macroscopic Maxwell 
equations (MMEs) under the explicit assumption that relevant 
time-averaged optical observables can be computed by 
averaging over a statistically uniform spatial distribution of all 
the particles throughout the volume V [1,4,5]. It is, therefore, 
imperative to examine whether polychromatic illumination 
(ie, averaging optical observables over a finite rather than 
infinitesimal spectral range) can be quantitatively equivalent to 
averaging over random particle coordinates (i.e., ensemble 
averaging). 


To the best of our knowledge this equivalence has never 
been demonstrated using controlled laboratory measurements. 
For practical reasons it may be more straightforward to give 
this demonstration using numerically exact computer solutions 
of the MMEs [7-9]. In fact, recent improvements in the 
superposition T-matrix method (STMM) [10,11] coupled with 
the current availability of efficient computer clusters make 
possible direct first-principle calculations of electromagnetic 
scattering by representative many-particle groups. Hence the 
main objective of this paper is to apply the modeling 
methodology developed in [12,13] to a quantitative analysis of 
whether the notion of electromagnetic scattering by a time- 
varying DRM can be applied to the case of polychromatic 
illumination of a fixed quasi-random particulate medium. 


2. MODELING METHODOLOGY 


For the purposes of our analysis, we will compare far-field 
scattering properties of two types of object. The first one 
(hereinafter Object 1) is a fixed N-particle configuration quasi- 
randomly populating an imaginary spherical volume having a 
radius R (Fig. 1). All particles have the same radius r and 
refractive index m; their coordinates are assigned using a 
random-number generator while making sure that the particle 
volumes do not overlap. The second type of object (hereinafter 
Object 2) is a spherical DRM whose scattering properties are 
modeled by taking the average of relevant optical observables 
over the equiprobable orientation distribution of Object 1. It 
has been demonstrated previously [5,12] that this type of 
averaging indeed yields numerical results that are essentially 
equivalent to those obtained by true ensemble averaging. 
Specifically, it has been shown that different realizations of 
Object 2 yield virtually indistinguishable quasi-monochromatic 
far-field observables. Another example of this equivalence will 
be discussed in Section 3. 


Fig. 1. A fixed N-particle configuration quasi-randomly populating an 
imaginary spherical volume. 
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where A is the wavelength, ni" = {6° ,@™°\ is the unit 
9 


vector in the incidence direction, n° = {6° ,g°@} is that in 
the scattering direction, p is the distance from the center of 
the object to the far-zone observation point, and {0,@} are the 
zenith and azimuth angles in the fixed (laboratory) spherical 
coordinate system centered at the object (Fig. 2). In what 
follows, we will assume that 6'"° = 0. 


We will consider three illumination scenarios. In the first 
scenario, Object 1 is illuminated by a quasi-monochromatic 


plane wave with a wavelength of Ag = 27/10 um. Note that the 


notion and physical causes of quasi-monochromaticity of light 
are discussed in [14]. The corresponding normalized scattering 
matrix is defined as 


FO) = COZ(O = 0, 9% = 0,0 = 0,9" =0;A), (2) 


where ch) is a normalization constant. 


Fig. 2. Far-zone scattering geometry. 


As is well known [4,5], the transformation of the time- 
averaged Stokes column vector of the incident quasi- 
monochromatic plane wave (“inc”) into the time-averaged 
Stokes column vector of the scattered spherical wavefront 
(“sca”) in the far zone of a fixed object is described by the phase 
matrix Z: 


In the second scenario, Object 1 is illuminated by a 
polychromatic plane beam consisting of quasi-monochromatic 
components with wavelengths uniformly distributed over the 
spectral range [0.95’9, 1.0529]. The 10% bandwidth is large 


enough to make the beam expressly polychromatic and yet 
small enough to be representative of many actual measurement 
conditions. The time-averaged Stokes column vectors of all the 
quasi-monochromatic components are assumed to be the same. 
Now the normalized scattering matrix is defined according to 


F2)@) = czas =@, os =0; ginc 7 0, gine =0;A))a, 
(3) 
where (---), denotes averaging over the entire spectral range 
in question. 
In the third scenario, Object 2 is illuminated by a quasi- 
monochromatic plane wave with the wavelength Apo. 


Furthermore, it is assumed that (i) the DRM represented by 
Object 2 is ergodic, and (ii) quasi-random oscillations of the 
incident electromagnetic field and the random changes in 
particle positions are statistically uncorrelated. According to 
Section 13.8 of [5], these assumptions imply that the 
corresponding time averaged normalized scattering matrix can 
be calculated as 
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Fig. 3. Elements of the matrices Fo) (thin black curves) and FA) (thick grey curves) for N = 200 and m = 1.32. 
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= cB)z(9s2 = @, ps; gine = 0, ig = gp; Ao; Py, 
(4) 


where (:::)y denotes ensemble averaging insofar as it is 
modeled by averaging over the uniform orientation distribution 
of Object 1 with respect to the laboratory coordinate system. 
Note that FO) is independent of o*“ because of taking the 
average over the uniform orientation distribution of Object 1. 

The normalization constants in Eqs. (2)-(4) are chosen such 
that the (1,1) element in each scattering matrix (i.e., the phase 
function) satisfies the normalization condition 


The scattering matrices F (6), F2)(@), and FO)(@) are 
computed using the numerically exact STMM solver of the 
MMEs described in [11]. In all computations, R is fixed at 5 ym 
and r is fixed at 0.4 um, which implies that the median size 
parameter of both types of particulate object is fixed at 
Xo = 2nR/Ao = 50 and that of the particles is fixed at 
x0 = 2ar/Ao = 4. Note that the STMM affords a quasi-analytical 
orientation-averaging procedure [10] that is highly efficient 
and accurate and yields FO)(@) results completely devoid of 
residual numerical noise (cf. [15]). The integration over the 
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Fig. 4. Panels (a), (b), and (c) depict the ratios FD(@)/FD(e) (thin black curves) and F2(0)/FO) (thick grey curves) for 
different realizations of Objects 1 and 2 (see text). Panels (d), (e), and (f) are the same as panels (a), (b), and (c), respectively, except 
that now the thin black curves depict the corresponding ratios F2(0)/FO (0) : 


spectral interval [0.9549, 1.0549] in Eq. (3) is performed by 


using the Gaussian quadrature formula with 100 division 
points. 


3. NUMERICAL RESULTS AND DISCUSSION 
Owing to the uniform orientation distribution of Object 1, the 


dimensionless scattering matrix F°)(6) has the following 
well-known symmetric structure [5,16]: 


F2) (0) FS) (®) F®) (0) F2) (0) 
(3 (3 (3 (3 
FO) -| OO) FPO) FPO) FPO) 6 
= FP) a Fe) Fe) FQ) 
FQ(@) FEe) -FP) Fee) 
with 
F8)(0) = F2)(0) = F(x) = F(x) = 0. (7) 


Our STMM computations for several realizations of Object 2 
have demonstrated that the elements populating the upper 
right and lower left 2x2 blocks of this matrix are negligibly 
small (in the absolute-value sense) compared to the other 
elements: 


FPO) FP) 0 0 
FO(@) ~ FO) FLO) 0 0 8 
7 EG) Z2@y{ (8) 
0 0) 6FX) FOO) 


0 0 6-F@(e) FE) 

One would expect this block-diagonal structure to hold exactly 
for a randomly oriented scattering object with a plane of 
symmetry [5,16]. Strictly speaking, a multi-particle object such 
as the one shown in Fig. 1 does not have a plane of symmetry. 
Yet by virtue of being composed of a large number of quasi- 
randomly positioned particles it probably possesses at least 
one quasi-plane of symmetry, which could explain Eq. (8). 

There is no a priori reason to expect the matrix FY) to 


have the same block-diagonal structure as well as satisfy Eq. 
(7). Figure 3 (computed for N = 200 and m = 1.32) shows 


indeed that the element Fe) does not vanish, the element 
FO) is distinctly different from the element FY), and 
there is little resemblance between the elements FO) and 
= FY). Furthermore, the elements Fe), Fe), 


Fo), and F{)(@) substantially deviate from zero at 


(000 qe = 100 
L E 4 4 
100 E 1 f 
E 4 50 - | 
10} | | 1 
E a| 0 =| J 
1E | | | 
E 4 -50 J ~ ig J 
0.1 | Fa3/F4, (%) 
E 4 | 
Qi ee -100 UO gy eh 
100 pow ————————— ——————— SSS 10 
50+ 4 sob 4 
oF ll ob 1 
OF FaalFrs (%) {OF —FaolF is (%) -Fy/F 14 (%) A 
199 44144 Pa ee J yop -..1.44 re pe J 
100 T T T T T T T T T T ry) © . oF 100 ——-——— >] 7] 7 100 ——-——— |] 7] 7 
50h 4 sok J 50h 4 
oF # of Opry WV oh 4 
“50r Fal Fas (%) 7 °F _FysfF yy (%) 7 OF FaalFar (%) 7 
oop. 41 sti. ii iit) -qo9l iii l peti iis J 499 1.1. La La J 
0 60 120 180 0 60 120 180 0 60 120 180 


Scattering angle (deg) 


Scattering angle (deg) 


Scattering angle (deg) 


Fig. 5. As in Fig. 3, but for the matrices F2)(6) (thin black curves) and FO) (thick grey curves). 


© =n. Most importantly, all Fe) curves in Fig. 3 reveal 
quasi-irregular large-amplitude oscillations which exemplify 
the famous speckle phenomenon typical of fixed objects 
illuminated by a collimated monochromatic or quasi- 
monochromatic beam [6,17]. 

In a stark contrast, the corresponding ensemble-averaged 


FA(@) curves satisfy Eq. (7) and—to a high numerical 
accuracy—Eq. (8), are completely devoid of speckles, and 
illustrate smooth scattering patterns typical of a true time- 
varying DRM [5,12]. 

Furthermore, a change in particle positions of Object 1 


changes the speckles in the FY (0) angular patterns 


completely, while leaving the corresponding matrix FS)(@) 


essentially intact. This is illustrated in Fig. 4(a) computed by 
switching the sign of the z-coordinates of the particles forming 
Object 1 in Fig. 3. In fact, the practical indistinguishability of all 


FI@) and FP O)/FPO) curves for this alternative 


realization of Object 2 from those for the original one (not 
shown here) illustrates once again that averaging over all 
orientation of a quasi-randomly generated Object 1 simulates 
very well the true ensemble averaging. 


Similarly dramatic changes in the FY) patterns occur 


when the number of particles increases from N = 200 to N = 
400 (Fig. 4(b)) or when the particle refractive index increases 
from m = 1.32 to m = 1.6 (Fig. 4(c)). Yet the corresponding 


FO)(0) curves remain quite smooth and speckle free. 


All these traits of the matrices FY) and F&)(0) show 


unequivocally that the latter is representative of a DRM 
whereas the former is not. 

Figure 5 reveals a dramatically different picture. It is seen 
that averaging over wavelengths has three pronounced effects. 
First, it drastically suppresses the speckles and yields much 
smoother scattering-matrix curves. Second, the structure of the 


matrix Fe) is now well described by Eqs. (7) and (8). 
Third, although there is residual small-amplitude “noise”, the 


overall angular profiles of the F2)(0) curves essentially 


coincide with the corresponding FO)(0) curves. These effects 


appear to be quite robust, as illustrated by comparing the top 
panels of Fig. 4 with their bottom counterparts. 

To explain the above STMM results qualitatively, let us 
assume that the spatial distribution of particles in Fig. 1 is so 
sparse that each particle is located in the far zones of all the 
other particles. In this case the exact integral-equation 
counterparts of the MMEs, called the Foldy equations, can be 
replaced by the algebraic far-field Foldy equations, as described 
in Chapter 6 of [5]. Then the scattered electromagnetic field in 
the far zone of the entire multi-particle group can be expressed 
as an infinite series of terms, each representing the 
contribution from an ordered n-particle sequence with n=1. 
The Stokes column vector (1) of the total scattered field can 
then be expressed as the sum of an infinite number of terms, 
each representing the result of interference of two spherical 
wavelets centered at the end particles of a pair of ordered 
particle sequences. 

Two such ordered particle sequences, with n = 4 and 3, are 
shown in Fig. 6. If the interference of the corresponding pair of 
spherical wavelets at the distant observation point is 
constructive (destructive) then it serves to increase (decrease) 


the total intensity scattered in the direction n°@. As a 
consequence, the pair of particle sequences contributes either a 
bright or a dark spot to the speckle pattern generated by the 
entire multi-particle group. 


Fig. 6. Interference origin of speckles generated by a sparse fixed 
multi-particle group. 


It is straightforward to show that the result of the 
interference of the wavelets generated by two ordered particle 
sequences depends largely on the complex exponential factor 
exp(iA), where i =(-1)'? and A is the phase difference 


between the corresponding light paths. In the case of the two 
light paths shown in Fig. 6, 


A= (22/A)\(r4 + R43 + R32 + Roy + nine Ry 
Ryy nize . Ry), (9) 


ry — R37’ 
where 7; is the distance from particle / to the observation point, 
Rj is the interparticle distance, and R; is the position vector of 


particle i with respect to the laboratory coordinate system. In 
general, random changes of particle coordinates in a DRM cause 
frequent temporal oscillations of the factor exp(iA). As a 


consequence, only the contributions from certain types of pairs 
of particle sequences survive time averaging and cause 
prominent optical phenomena typical of DRMs such as the 
forward-scattering interference, the diffuse background, and 
coherent backscattering [5,12]. All three phenomena are best 


seen in the resulting phase function F(0) (the upper left 


panel of Fig. 3) in the form of the strong forward-scattering 
peak, featureless side scattering, and a narrow backscattering 
maximum. 

In a fixed multi-particle configuration such as that shown in 
Fig. 1, the coordinates of the particles do not change. However, 
it is obvious from Eq. (9) that another efficient way of 
randomizing the factor exp(iA) is to scan a range of 


wavelengths. This explains qualitatively why the wavelength- 
averaged matrix F2)(@) in Figs. 4(d)-4(f) and 5 exhibits the 
same overall traits as those typical of the ensemble-averaged 


matrix F°)(6) computed for a true time-varying DRM. 
The residual small-amplitude noise in the FO) curves 


and the complete lack thereof in the F)(@) curves could be 
explained by the smaller average range of A/2z values in the 


case of averaging over wavelengths for a relatively small 
particulate volume. Indeed, it is a large number of oscillations 
of the factor exp(iA) that serves to reduce the residual noise. If 


so, this noise can be expected to weaken and _ ultimately 
disappear with increasing Xo and N. 


4. CONCLUDING REMARKS 


In view of Eq. (9) and the accompanying discussion, it may not 
be exceedingly surprising that a fixed, quasi-random particulate 
medium illuminated by a polychromatic parallel beam exhibits 
qualitatively the basic scattering traits of a true DRM. What is 
surprising is the numerical near identity of the respective 
scattering patterns despite the large values of the packing 
density (10% for N = 200 and 20% for N = 400). The numerical 
closeness of the FC)(@) and F°)(0) curves in Figs. 4(d)-4(f) 
and Fig. 5 cannot be explained by simple qualitative arguments 
based on the far-field Foldy equations and could only be 
revealed via direct, numerically exact solutions of the MMEs. 
Yet this result appears to be fundamentally important since it 
allows for the generalization of the notion of electromagnetic 
scattering by a DRM to encompass not only time-varying 
particulate media but also fixed quasi-random multi-particle 
configurations provided that the latter are illuminated by 
polychromatic light. 

The suppression of speckles in Figs. 4(d)-4(f) and Fig. 5 
appears to give credence to the widespread use of the radiative 
transfer theory and its ad hoc modifications to model “diffuse 
multiple scattering” by fixed particulate layers illuminated by 
polychromatic beams. One should remember however that this 


theory is explicitly based on the far-field Foldy equations [4,5] 
and as such may be inapplicable to particulate samples with 
packing densities exceeding a few percent [18-20]. 

The 10% bandwidth of the polychromatic beam used in our 
simulations can be larger than that in many actual experiments. 
However, the typical size parameter of real particulate samples 
is also much greater than 50, which implies that the resulting 
range of variability of A in Eq. (9) should be sufficiently large 
for our main conclusions to hold. 

Finally we note that polychromatic illumination is not the 
only factor potentially causing the suppression of the observed 
speckle pattern. Among other factors are the use of detectors of 
light with poor angular resolution and illumination by 
uncollimated incident beams [6] (e.g., Gaussian beams [21]). It 
should be quite instructive to extend our work by modeling the 
effects of these factors—as well as their combinations with 
spectral averaging—on the speckle patterns generated by fixed 
particulate samples. Yet another factor can be a macroscopic 
movement of a rigid particulate sample as a whole relative to 
the source of light and/or the detector during the 
measurement. This movement also can result in dramatic 
quasi-random changes of the exp(iA) factors and cause the 


suppression of the resulting speckle pattern, thereby 
revealing—at least in a qualitative sense—the typical diffuse 
and coherent-backscattering regimes [22-24]. This last factor 
can be especially important in remote-sensing applications 
(see, e.g., [25-27]). 
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